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摘 要 : 研究 融雪 末期 雪线 高 度 变化 有 助 于 预测 冰雪 系统 未 来 的 变化 趋势 .理解 区 域 和 全 球 气候 变化 。 基 于 
Google Earth Engine (GEE ) 平 台 和 Landsat 卫星 数据 ,发 展 了 区 域 雪线 高 度 提 取 模 型 ,提取 了 天 山 4 个 子 流域 1991 一 


2021 年 的 融雪 末期 雪线 高 度 ,并 分 析 了 雪线 高 度 的 变化 特征 及 其 与 气象 因素 的 关系 。 结 果 表 明 : (1) Landsat 提 取 
的 融雪 末期 雪线 高 度 与 Sentinel-2 提 取 的 融雪 未 期 “最 小 化 " 积 雪 范围 具有 较 高 的 一 致 性 ,其 总 体 精度 为 91.6%， 


Kappa 系数 达 0.9 以 上 ,利用 该 模型 可 准确 获取 融雪 末期 


的 区 域 雪 线 高 度 ,(2) 研究 区 近 30 a 融雪 末期 雪线 高 度 呈 明 


显 的 上 升 趋势 ,上 升 速率 介 于 2.7~6.4 ma 之 间 ; 其 中 , 玛 纳 斯 河流 域 雪线 高 度 上 升 速度 最 快 , 阿 克 牙 孜 河流 域 雪线 


高 度 上 升 速度 最 慢 。(3) 夏季 气温 是 影响 研究 区 融雪 末 
BUB 


期 雪线 高 度 变化 的 主要 因素 (P< 0.05) ,降水 对 其 影响 相对 


关键 词 : SARE; Landsat; 遥感 监测 ; 气候 变化 ; KU 


雪线 作为 地 球 陆地 表面 积 雪 覆盖 区 与 非 积 雪 
履 盖 区 之 间 的 分 界线 ,其 高 度 变 化 信息 是 冰雪 水 
资源 变化 的 直观 反映 。 在 中 纬度 地 区 ,季节 变化 能 
引起 雪线 位 置 的 升降 ,这 种 临时 界线 也 称 为 季节 雪 
线 ( 或 瞬时 雪线 )2 5 。 而 在 冰川 学 中 雪线 特 指 消融 
末期 积 雪 存 在 的 下 限 外 ,等 同 于 融雪 末期 雪线 “7” 。 
融雪 末期 雪线 相 比 季节 雪线 更 具有 气候 意义 ,常用 
于 指示 冰川 物质 平衡 < … ,直接 反映 着 冰川 的 进退 ; 
研究 融雪 末期 雪线 高 度 变 化 有 助 于 预测 冰川 和 非 
季节 性 积 雪 未 来 的 变化 趋势 .理解 区 域 和 全 球 气候 
变化 路 。 因 此 ,准确 获取 融雪 末期 雪线 高 度 信息 并 
实现 其 动态 监测 ,不 仅 可 为 研究 冰冻 圈 变 化 及 其 对 
气候 变化 的 响应 提供 关键 信息 ,对 寒 区 水 资源 管理 
及 其 可 持续 利用 也 具有 重要 意义 。 

遥感 技术 在 冰冻 圈 要 素 的 时 空 变化 研究 中 发 
挥 着 重要 作用 ,已 成 为 冰冻 圈 研 究 的 重要 技术 手 
段 ”“。MODIS 表 感 数据 因 高 时 间 分 辨 率 的 特点 ， 


融雪 末期 雪线 高 度 遥 感 提取 模型 ,提高 了 对 亚洲 高 
山区 雪线 高 度 的 时 空 分 布 规律 的 宏观 认识 ;但 MO- 
DIS 数据 因 受 空间 分 辩 率 的 限制 , 仅 适 用 于 大 尺度 
融雪 末期 雪线 高 度 的 宏观 遥感 监测 。Landsat 系列 
卫星 提供 了 长 时 序 和 较 高 分 辨 紊 的 各 感 数据 (近年 
发 展 的 Sentinel-2 卫星 时 空 分 辩 率 更 高 ) ,然而 传统 
计算 机 难以 满足 处 理 海量 遥感 数据 所 需 的 存储 与 计 
算 资源 ,使 得 这 些 高 分 辨 率 的 遥感 大 数据 尚未 被 充 
分 挖掘 和 应 用 。 现 有 基于 Landsat 系列 遥感 数据 的 
雪线 高 度 监测 研究 仅 限于 小 范围 局 地 区 域 “”*。 
近年 来 ,遥感 云 计算 平台 Google Earth Engine (GEE) 
的 出 现 为 遥感 大 数据 挖掘 提供 了 条 件 , 并 改变 了 传 
统 遥 感 数据 处 理 和 分 析 的 模式 ,使 得 全 球 范 围 长 时 
间 序 列 的 分 析 和 应 用 成 为 可 能 ;GEE 强大 的 存储 
与 云 计算 能 力 在 积 雪 脖 感 领域 也 呈现 出 特有 的 优 
Bs 2028) 。 此 外 ,如 何 降低 光学 遥感 中 云 诞 蔽 的 影 
响 ,一 直 是 积 雪 遥感 监测 中 面临 的 重要 问题 ,也 是 


在 大 范围 季节 (或 瞬时 ) 雪 线 高 度 监测 中 应 用 |) 
iz 79579, Tang 等 发展 了 基于 MODIS 积 雪 产 品 的 
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有 效 获取 雪线 高 度 信息 的 关键 。 因 此 ,综合 高 分 辩 
率 遥 感 大 数据 的 优势 ,降低 云 遮 蔽 的 影响 ,发 展 遥 


点 实验 室 开 放 基 金 项 目 (SKLCS-OP-2020-08 ) ;湖南 省 教育 厅 科研 项 目 


(20B227) 


作者 简介 : 王 靖 文 (1998-) , 男 ,硕士 研究 生 ,主要 从 事 积 雪 遥 感 方 


的 研究 . E-mail: wangjingwen@mail.hnust.edu.cn 


通讯 作者 : 唐 志 光 . E-mail: tangzhgl1@hnust.edu.cn 


1385 一 1397 页 http: //azr.xjegi.com 


202212.00143v1 


chinaXiv 


1386 Tom 


区 研 R 


ChinaXiv 合 作 期 刊 


39 卷 


感 云 计算 平台 支持 下 的 融雪 末期 雪线 高 度 提 取 模 
型 ,是 准确 监测 大 范围 .长 时 序 融 雪 末 期 雪线 高 度 
变化 的 重要 途径 。 

天 山 被 誉 为 "中 亚 水 塔 ”, 横 百 于 欧 亚 大 陆 腹 
地 ,由 一 系列 高 大 山地 . 山 间 盆地 和 和 谷地 组 成 ,是 世 
界 上 距离 海洋 最 远 .现代 冰川 发 育 最 为 广泛 的 山系 
之 一 下 。 天 山 丰 富 的 冰川 . 积 雪 是 中 亚 干旱 区 的 重 
要 水 资源 ,对 气候 变化 非常 敏感 。 在 过 去 的 半 个 多 
世纪 (1960 一 2016 年 ) ,天 山 的 年 均 气 温 呈 显著 增加 
趋势 , 增 速 为 0.32 C- (10a), 比 全 球 气温 增加 速率 
高 出 0.19 C- (103) 99, fee AS BE DEL RL KS 
变化 直接 影响 河川 径流 过 程 与 水 资源 量 的 改变 , 导 
致 中 亚 地 区 水 资源 的 供需 矛盾 加 剧 ,影响 中 亚 区 域 
国家 之 间 关 系 以 及 “一 带 一 路 ”倡议 ”。 为 此 ,开展 
天 山 雪线 高 度 变 化 的 遥感 监测 及 其 气候 响应 研究 ， 
可 为 研究 气候 变化 对 天 山 山区 的 水 资源 影响 提供 
重要 信息 ,对 科学 管理 水 资源 .服务 于 "一 带 一 路 ” 
倡议 具有 重要 意义 。 

本 文选 取 天 山地 区 的 4 个 山区 流域 作为 研究 对 
象 , 以 Landsat 系列 (TM ETM+ , OLD 卫星 数据 为 主 
要 数据 源 , 在 遥感 云 计算 平台 GEE 支持 下 ,发 展 了 
融雪 末期 雪线 高 度 提取 模型 ,提取 了 4 个 流域 
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(c) 木 扎 提 河 流域 
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1991—2021 年 的 融雪 末期 雪线 高 度数 据 ,并 分 析 了 
人 研究 区 融雪 末期 雪线 高 度 的 变化 特征 及 其 与 气候 
变化 的 关系 。 


1 研究 区 概况 及 数据 


1.1 研究 区 概况 

天 山 山 脉 横 跨 中 国 、 哈 萨克斯 坦 .吉尔吉斯 坦 
和 乌兹别克 斯 坦 等 国家 。 天 山地 形 复杂 ,高 山区 降 
水 量 大 且 气 温 低 ,冰川 发 育 广泛 ,其 冰雪 变化 关系 
着 诸多 流域 的 水 资源 环境 与 水 文 过 程 。 本 文选 择 
琼 乌 散 库 什 河 流域 . 木 扎 提 河流 域 . 阿 克 牙 孜 河流 
域 和 玛 纳 斯 河流 域 作为 研究 的 典型 流域 (图 1), 它 
们 均 发 源 于 冰雪 覆盖 区 ,冰雪 融 水 在 河川 径流 中 占 
有 重要 比例 。 琼 乌 散 库 什 河和 木 扎 提 河 发 源 于 天 
山南 坡 , 最 终 均 汇 人 塔里木 河流 域 , 琼 乌 散 库 什 河 
流域 海拔 4000 m 以 上 ,主要 是 冰雪 和 覆盖, 约 占 整 个 
流域 面积 的 35% , 木 扎 提 河 流域 内 发 育 有 大 量 的 现 
代 冰 川 , 约 占 整 个 流域 面积 的 48%; 阿 克 牙 孜 河流 
域 约 20% 集 中 在 海拔 高 程 4000 m 以 上 的 高 山 带 ; 玛 
纳 斯 河 发 源 于 天 山北 坡 , 流 入 准噶尔 盆地 腹地 , 流 
域内 海拔 4000 m 以 上 常年 被 永川 和 非 季 节 性 积 雪 
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(e) 玛 纳 斯 河流 域 


图 1 研究 区 概况 图 
Fig. 1 Overview of the study area 
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TEE o , 约 占 整个 流域 面积 的 15%。 
1.2 数据 来 源 

研究 采用 了 Landsat TM/ETM+/OLI、 Sentinl- 2 
MSI,STRM DEM 和 ERA5-LAND 再 分 析 数 据 集 ,这 
些 数据 均 通 过 CEE 云 计算 平台 在 线 访问 和 处 理 。 
1.2.1 Landsat 数 据 采用 Landsat 系 列 数据 中 的 陆 
表 反 射 率 (Land Surface Reflectance ) 数 据 ,该 数据 已 
经 过 辐射 定 标 和 大 气 校正 ,消除 了 大 气 散射 .吸收 、 
反射 引起 的 误差 ,空间 分 辨 率 为 30 m, 时 间 分 辨 率 
为 16 d。 该 数据 提供 了 由 CFMask 去 云 算法 生成 的 
质量 评估 (QA) 波 段 , 包 含 了 对 云 和 云 阴影 的 标识 。 
图 2 统计 了 各 研究 区 逐年 融雪 期 (6 一 9 月) 共享 的 
Landsat 影像 数量 ,以 及 研究 区 内 云 量 小 于 20% 的 影 
像 频次 ;利用 Landsat 数 据 可 以 实现 研究 区 1991 年 
至 今 的 融雪 末期 雪线 高 度 监测 。 
1.2.2 Sentinel-2 MSI 数 据 欧洲 航天 局 (ESA) 发 射 
的 Senitnel-2 卫星 由 Sentinel-2A 和 Sentinel-2B 组 
成 ,分 别 于 2015 年 ,2017 年 发 射 ,双星 重 访 周期 为 5 
d, 均 携带 多 光谱 成 像 仪 (MSI) ,共有 13 个 通道 ,空间 
分 辨 率 有 10 m,20 m 和 60 m。 本 文采 用 Sentinel-2 
MSI 数 据 的 Level-2A 产 品 对 Landsat 数据 的 雪线 高 
度 提取 结果 进行 验证 ;该 产品 为 经 过 辐射 定 标 和 大 
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气 校正 的 地 表 反 射 率 产 品 , 选 
辩 率 为 10 mo 

1.2.3 SRTM DEM 数据 ”用 于 提取 雪线 高 度 信息 的 
DEM 数据 为 SRTM DEM 数据 (NASA 3.0 版 )。 该 数 
据 通过 ASTER GCDEM2 ,GMTED2010 和 NED 数据 对 
原始 SRTM 数 据 的 缺失 值 进行 填补 ,空间 分 辨 率 为 
30 m。 该 数据 的 全 球 垂直 精度 约 为 16 m, 在 中 国 地 
形 平坦 地 区 ,绝对 垂直 精度 约 为 14 m, 地 形 起 伏 较 
大 山区 绝对 垂直 精度 约 为 26 mo 

1.2.4 ERAS-LAND 再 分 析 数 据 集 ERAS- LAND 
再 分 析 数 据 集 来 源 于 欧洲 中 期 天 气 预 报 中 心 (Euro- 
pean Centre For Medium- Range Weather Forecasts, EC- 
MWF) ,该 数据 集 提 供 了 1979 年 1 月 至 今 的 全 球 降 
水 .气温 等 变量 数据 。 本 文采 用 该 数据 集中 的 月 平 


用 数据 波段 的 空间 分 


均 气 温和 降水 数据 ,空间 分 辩 率 为 10 km, 以 探究 研 
究 区 融雪 末期 雪线 高 度 的 变化 与 气候 因子 之 间 的 
响应 关系 。 


2 研究 方法 


选取 接近 融雪 末期 的 多 时 相 Landsat TM/ETM+/ 
0LI 影 像 为 数据 源 , 构 建 雪 线 高 度 提 取 模 型 的 技术 


路 线 如 图 3 所 示 。 具 体 的 方法 主要 包括 积 雪 遥感 识 


= Landsat5 = Landsat7 = Landsat 8 云 量 小 于 20% 


80 「 (a) 琼 乌 散 库 什 河流 域 


影像 数量 
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(c) 阿 克 牙 孜 河流 域 


影像 数量 
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50 - (b) 木 扎 提 河 流域 
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(d) 玛 纳 斯 河流 域 
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图 2 研究 区 域 融 雪 期 覆盖 的 影像 数量 


Fig.2 Number of image covered by the study area at the melting Season 
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接近 融雪 未 期 的 多 时 相 
Landsat TM/ETM+/OLI 


多 时 相 积 雪 分 类 图 像 
( 积 轨 、 云 、 神 地 ) 


区 域 雪线 高 度 
去 云 、 时 序 融 合 


融雪 末期 的 “最 小 化 " 积 雪 分 布 图 


--[ 精度 评价 


30m DEM 


融雪 末期 雪线 高 度 提 取 


1991 一 2021 年 融雪 末期 雪线 高 度 变化 


图 3 总 体 技术 路 线 


Fig.3 Overall technical route 


别 .融雪 末期 "最 小 化 " 积 雪 范围 提取 、 区 域 雪 线 高 
度 提 取 和 精度 验证 。 
2.4 积 雪 遥感 识别 
采用 SNOWMAP 算 法 ”提取 积 雪 范 围 。 该 算 

法 基于 积 雪 在 可 见 光波 段 , 具 有 较 高 的 反射 率 ,在 
短 红 外 波段 具有 较 强 的 吸收 特征 ,建立 归 一 化 积 雪 

EJ (Normalized Difference Snow Index, NDSI) ,设置 
合适 的 国 值 ,用 以 确定 积 雪 覆盖 范围 。NDSI 计 算 公 
式 为 : 


NDST= P Green 一 PSswm (1) 
P Gees + Pswik 


式 中 :pe 和 pswn 分 别 代 表 Landsat 数据 中 的 绿 光波 
段 和 短波 红外 的 反射 率 值 。 由 于 水 体 和 积 雪 具有 
相似 的 波段 反射 特性 ,为 了 避免 水 体 .阴影 被 误 识 
别 为 积 雪 ,增加 1 个 判别 条 件 即 近 红外 波段 反射 率 
值 (pur)>0.11。 相 关 研 究 证 明了 NDSI [IS] f& y 0.29 
更 适合 于 青藏 高 原 地 区 的 积 雪 制图 六.。 因 此 ,本 研 
究 将 NDSI> 0.29 & pw>0.11 的 像 元 识别 为 积 雪 像 
元 ,不 满足 上 述 条 件 的 归 类 为 非 雪 像 元 ,最终 将 识 
别 得 到 的 积 雪 , 非 雪 与 OA 波段 识别 得 到 的 云 合成 
得 到 积 雪 分 类 图 ( 雪 、 云 非 雪 )。 
2.2 融雪 末期 “最 小 化 ? 积 雪 范 围 提 取 

通过 瞬时 雪线 高 度 去 云 和 时 序 融 合 ,以 减少 云 
的 影响 并 得 到 最 能 代表 融雪 末期 的 “最 小 化 " 积 雪 


2.2.1 瞬时 雪线 高 度 去 云 ” 对 区 域内 云 和 雪 的 比例 
进行 判断 ,车 云 和 雪 面 积 比例 <1/2, 则 提取 该 区 域 
的 瞬时 雪线 高 度 以 进行 去 云 处 理 。 瞬 时 雪线 去 云 
策略 :如 果 云 像 元 对 应 的 地 面 海拔 高 度 大 于 瞬时 区 
域 雪 线 高 度 , 则 归 类 为 积 雪 ;反之 , 则 归 类 为 陆地 。 
因 云 覆盖 的 比例 过 大 ,容易 出 现 积 雪 像 元 极 少 或 者 
积 雪 履 盖 区 周围 均 为 云 覆 盖 ,导致 提取 的 瞬时 雪线 
高 度 存 在 较 大 偏差 或 者 无 法 提取 。 为 了 保证 提取 
的 瞬时 区 域 雪线 高 度 的 有 效 性 ,本 研究 设 定 云 和 雪 
面积 比例 >1/2 时 , 则 直接 参与 时 序 融 合 处 理 。 
2.2.2 WERA ”对 接近 融雪 末期 的 多 时 相 积 雪 分 
类 图 像 ( 包 含 瞬 时 雪线 高 度 去 云 后 图 像 ) 进 行 时 序 
融合 ,以 完成 积 雪 范围 的 最 小 化 合成 ,最 终 得 到 能 
代表 融雪 末期 的 “最 小 化 ” 积 雪 范 围 图 。 时 序 融 合 
策略 :对 每 一 像 元 ,如 果 被 任 一 时 相识 别 为 裸 地 , 则 
归 类 为 裸 地 ;如 果 50% 以 上 的 像 元 被 识别 为 积 雪 ， 
而 其 他 时 相 都 为 云 , 则 归 类 为 积 雪 ;如 果 50% 以 上 
的 像 元 被 识别 为 云 ,而 其 他 时 相 都 为 雪 , 则 归 类 为 
云 。 融 合 后 得 到 的 代表 融雪 末期 的 “最 小 化 " 积 雪 
范围 图 ,用 于 融雪 末期 雪线 高 度 的 提取 。 
2.3 区 域 雪线 高 度 (RSLA) 提 取 

区 域 雪 线 高 度 法 的 原理 为 :在 目标 区 域内 找 出 
一 个 高 程 , 使 得 该 高 程 之 上 的 非 雪 像 元 数量 与 该 高 
程 之 下 的 积 雪 像 元 数量 之 和 最 小 ,那么 这 一 高 程 就 
代表 该 区 域内 的 雪线 高 度 (图 4) ,公式 表示 为 : 

sum = Cy, + Ca (2) 


式 中 :sum 为 高 程 之 下 积 雪 像 元 数量 (CG,) 与 高 程 之 
上 非 雪 像 元 数量 (C6) 之 和 ;RSLA 为 区 域 雪线 高 
度 ;elevii 为 区 域内 最 低 高 程 ;elev 为 区 域内 最 高 
高 程 。 

具体 过 程 为 :从 区 域 最 低 高 程 开始 ,计算 目标 
高 程 之 下 的 积 雪 像 元 数量 与 该 高 程 之 上 的 非 雪 像 
元 数量 之 和 ,然后 每 隔 1 mm 高 程 进行 1 次 该 计算 , 直 
到 区 域内 最 高 高 程 。 然 后 利用 最 小 值 函 数 (Min ) 来 
检索 数量 和 中 的 最 小 值 sumas , sumas PT XT D. AY tes 
程 即 为 目标 区 域 的 雪线 高 度 。 该 方法 可 以 允许 部 
分 云 履 盖 的 存在 ,是 实现 目标 区 域 雪线 高 度 的 最 优 
估算 ,因而 也 适用 于 积 雪 图 像 的 去 云 处 理 。 
2.4 雪线 高 度 的 精度 评估 

采用 空间 分 辨 率 更 高 . 重 访 周 期 更 短 的 Senti- 
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81°10'E 81°40'E 100 
(a) 阿 克 牙 孜 河流 域 9 高程 之 下 的 积 雪 像 元 数 
a 高 程 之 上 的 陆地 像 元 数 
Z 
F+ g É 
S e 
FERS 
z "$^ jz 
EI 一 区 域 雪线 高 度 (RSLA)| 和 
vicum mimozmez]|* | 
81?10'E 81°40'E min RSLA max 
图 4 区 域 雪线 高 度 提取 示意 图 
Fig.4 Sketch map of the extraction of regional snowline altitude 
nel-2 影 像 ,通过 提取 融雪 末期 “最 小 化 " 积 雪 范围 ， Ree= LA x100 (6) 
十 


来 验证 Landsat 数据 提取 融雪 来 期 区 域 雪 线 高 度 的 
精度 。 具 体验 证 方法 :选取 2019 一 2021 年 研究 区 接 
近 融 雪 末 期 且 云 量 低 于 10% 的 多 时 相 Sentinel-2 遥 
感 影像 ( 表 1) ,采用 目 视 解 译 法 ,对 每 一 影像 进行 积 
雪 范 围 制图 ;然后 将 接近 融雪 末期 的 多 时 相 积 雪 范 
围 图 进行 时 序 融 合 , 得 到 Sentinel-2 提 取 的 融雪 末 
期 “最 小 化 " 积 雪 范 围 (作为 融雪 末期 积 雪 范围 “ 真 
值 ”)。 最 后 ,基于 Sentinel-2 提 取 的 融雪 末期 “最 小 
化 ”" 积 雪 范 围 图 ,在 雪线 位 置 附近 进行 验证 样本 的 
随机 采样 ,并 结合 DEM 数据 和 Landsat 提 取 的 融雪 
末期 区 域 雪 线 高 度 值 ,计算 融雪 末期 雪线 高 度 提取 
模型 的 总 体 精度 (0A) .精确 率 (Pre) 和 召回 率 
(Rec) ,计算 公式 如 下 : 

OA= 


A+D 
A+B+C+D 


x 100 (5) 


x 100 (4) 


A 
Pre = 
we A+B 
R1 研究 区 域 Sentinel-2 影像 覆盖 数量 
Tab.1 Number of Sentinel-2 image covered by 
the study area 
TE ACE 
琼 乌 散 库 什 河流 域 26 


分 布 时 间 / 年 -月 -日 
2019-06-01 一 2019-09-30 
2020-06-01 一 2020-09-30 
2021-06-01 一 2021-09-30 
2019-06-01 一 2019-09-30 
2020-06-01 一 2020-09-30 
2021-06-01 一 2021-09-30 
2019-06-01 一 2019-09-30 
2020-06-01 一 2020-09-30 
2021-06-01 一 2021-09-30 
2019-06-01 一 2019-09-30 
2020-06-01 一 2020-09-30 
2021-06-01 一 2021-09-30 


木 扎 提 河 流域 24 


阿 克 牙 孜 河流 域 26 


玛 纳 斯 河流 域 24 


式 中 :4 为 高 于 Landsat 融 雪 末 期 雪线 高 度 且 在 Sen- 
tinel-2“ 最 小 化 " 积 雪 范围 中 识别 为 积 雪 的 点 数 ;B 
为 高 于 Landsat 融雪 末期 雪线 高 度 且 在 Sentinel-2 
“最 小 化 " 积 雪 范围 中 识别 为 非 雪 的 点 数 ;C 为 低 于 
Landsat 融雪 末期 雪线 高 度量 在 Sentinel-2“ 最 小 化 ” 
积 雪 范围 中 识别 为 积 雪 的 点 数 ;D 为 低 于 Landsat 融 
雪 示 期 雪线 高 度 且 在 Sentinel-2“ 最 小 化 ” 积 雪 范 围 
中 识别 为 非 雪 的 点 数 。Pre>Rec 时 ,代表 雪线 高 度 
存在 高 估 ;Pre<Rec 时 ,代表 雪线 高 度 存在 低估 。 


3 BRAD 


3.1 融雪 末期 雪线 高 度 精度 评估 结果 分 析 

通过 建立 混 消 和 矩阵 ,计算 研究 区 域 融雪 末期 雪 
线 高 度 提取 结果 的 Pre、Rec、OA 和 Kappa 系 数 ,对 提 
取 结 果 的 精度 进行 定量 评价 ( 表 2)。 结 果 表 明 :4 个 
流域 融雪 未 期 雪线 高 度 提 取 结 果 具 有 较 高 的 精度 ， 
OA 介 于 87.5%~94.4% 之 间 ,Kappa 系数 介 于 0.86~ 
0.93 之 间 ;平均 OA 和 Kappa 系数 分 别 为 91.6% 和 
0.9。 然 而 ,4 个 流域 Pre 均 略 低 于 Rec, 说 明 相 比 时 
间 分 辩 率 更 高 的 Sentinel-2 融 雪 末 期 “最 小 化 ” 积 雪 
范围 ,Landsat 提 取 的 融雪 末期 雪线 高 度 略 有 低估 。 
以 阿 克 牙 孜 河流 域 为 例 ,0A 为 92.7% ,说 明 Landsat 
提取 的 融雪 末期 雪线 高 度 与 Sentinel-2 提 取 的 融雪 
末期 “最 小 化 ” 积 雪 范围 之 间 具 有 较 高 的 一 致 性 ; 精 
确 率 Pre 为 87.06% ,说 明 在 Landsat 提 取 的 融雪 末期 
雪线 高 度 之 上 日 在 Sentinel-2 提 取 的 “最 小 化 ” 积 雪 
范围 中 被 识别 为 非 积 雪 的 概率 为 12.94% ;而 该 流域 
Rec 为 93.97% ,说 明 在 Landsat 提取 雪线 高 度 之 下 且 
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表 2 Landsat 融 雪 末 期 雪线 高 度 提取 结果 精度 评价 
Tab. 2 Accuracy evaluation of snowline altitude extraction results at the end of melting season 
Sentinel-2 提 取 的 验证 样本 
didi 积 雪 /个 数 非 雪 /个 数 ani 
琼 乌 散 库 什 河流 域 高 于 Landsat 提 取 雪 线 高 度 291 76 Pre=79.29% 
IRF Landsat ERER HE 49 584 Rec=85.59% 
OA=87.5% 
Kappa=0.86 
木 扎 提 河 流域 高 于 Landsat 提 取 雪 线 高 度 391 68 Pre=85.19% 
低 于 Landsat 提 取 雪 线 高 度 16 525 Hee 
OA=91.6% 
Kappa=0.90 
阿 克 牙 孜 河流 域 高 于 Landsat 提 取 雪 线 高 度 343 51 Pre=87.06% 
低 于 Landsat 提 取 雪 线 高 度 22 584 Rec=93.97% 
OA=92.7% 
Kappa=0.92 
玛 纳 斯 河流 域 高 于 Landsat 提 取 雪 线 高 度 406 34 Pre=92.27% 
低 于 Landsat 提 取 雪 线 高 度 22 517 Rec=94.86% 
OA=94.4% 
Kappa=0.93 


在 Sentinel-2 提 取 的 “最 小 化 " 积 雪 范 围 中 被 识别 为 
积 雪 的 概率 仅 为 6.03%。 
3.2 融雪 期 雪线 高 度 时 空 变 化 特征 

获取 各 流域 6 一 9 月 以 16 d 为 间隔 的 多 年 雪线 
高 度数 据 (图 5)。 雪 线 高 度 在 6 一 9 月 呈现 出 强烈 的 
季节 性 变化 。 从 6 月 开始 , 随 着 积 雪 的 快速 融化 ,人 研 
究 区 积 雪 覆 盖 率 逐渐 降低 .雪线 高 度 逐 渐 升 高 , 瞬 
时 雪线 高 度 的 最 大 值 主要 出 现在 8 月 ,而 8 月 以 后 
雪线 高 度 又 逐渐 降低 ; 且 研 究 区 瞬时 雪线 高 度 体现 
出 较 强 的 年 际 波动 (图 5)。 以 2020 年 阿 克 牙 孜 河流 
域 的 瞬时 雪线 高 度 变化 为 例 ( 图 6) ,雪线 高 度 在 6 


域 近 30 a 融雪 末期 雪线 高 度 上 升 约 135 m; AFL 
提 河 流域 多 年 平均 融雪 末期 雪线 高 度 为 4207 m, 趋 
势 斜率 为 3.1 ma”, UE 30 a 融雪 末期 雪线 高 度 上 升 
2593 m; 玛 纳 斯 河流 域 融雪 末期 雪线 高 度 总 体 上 与 
木 扎 提 河 流域 接近 ,该 流域 多 年 平均 融雪 末期 雪线 
高 度 为 4212 m, 趋势 斜 率 最 高 ,为 6.4m*a', 近 30 a 
融雪 末期 雪线 高 度 上 升 约 192 m; 阿 克 牙 孜 河流 域 
多 年 平均 融雪 末期 雪线 高 度 为 3955 m ,趋势 斜率 在 
4 个 流域 中 最 低 ,为 2.8 ma”, JE 30 a 融雪 未 期 雪线 
高 度 上 升 约 84 m。 此 外 ,通过 计算 研究 区 融雪 末期 
雪线 高 度 的 标准 差 (o ) ,来 反映 融雪 末期 雪线 高 度 


8 月 逐渐 升 高 ,从 6 月 的 3775 m 上升 到 8 月 的 3936 
m, 而 9 月 雪线 高 度 又 下 降 到 3458 m; 经 模型 提取 的 
融雪 末期 雪线 高 度 为 3967 m, 高 于 每 一 瞬时 雪线 
高 度 。 
3.3 融雪 末期 雪线 高 度 年 际 变化 特征 

1991—2021 年 研究 区 域 融 雪 末 期 雪线 高 度 总 
体 上 呈 明 显 的 上 升 趋势 (图 7)。 琼 乌 散 库 什 河流 域 
和 玛 纳 斯 河流 域 融 雪 末 期 雪线 高 度 最 高 值 均 出 现 
在 2018 年 ,最 低 值 分 别 出 现 在 2003 年 和 1994 年 ; 木 
扎 提 河流 域 和 阿 克 牙 孜 河 流域 融雪 末期 雪线 高 度 
最 高 值 最低 值 分 别 出 现 在 2016 年 .2015 年 和 1993 
年 .1998 年 ,4 个 流域 近 30 a 融雪 末期 雪线 高 度 振幅 
均 不 超过 400 m。 琼 乌 散 库 什 河流 域 多 年 平均 融雪 
末期 雪线 高 度 为 4504 m ,趋势 斜率 为 4.5 mea! AU 


的 年 际 波动 幅度 (图 7)。 结 果 表明 : 木 扎 提 河 流域 
融雪 末期 雪线 高 度 标准 差 最 小 (o=59 m) ,说 明 该 流 
域 的 融雪 末期 雪线 高 度 年 际 波动 最 小 , 芍 乌 散 库 什 
河流 域 融雪 末期 雪线 高 度 标准 差 最 大 (o=94 m) ,说 
明 该 流域 对 应 的 融雪 末期 雪线 高 度 年 际 波动 最 大 。 
3.4 融雪 末期 雪线 高 度 变 化 与 气温 、 降 水 之 间 的 
关系 

雪线 高 度 与 气候 之 间 的 关系 非常 密切 。 本 文 
通过 计算 年 平均 和 夏季 平均 气温 .年 总 和 夏季 总 降 
水 量 与 融雪 末期 雪线 高 度 之 间 的 相关 关系 ,分析 融 
雪 末 期 雪线 高 度 与 气候 变化 之 间 的 响应 关系 ( 表 
3) ,其 中 ,年 降水 量 和 年 气温 指 的 是 由 上 一 年 10 月 
到 当年 9 月 的 总 降水 量 和 平均 气温 。 琢 乌 散 库 什 河 
流域 年 降水 量 和 和 夏季 降水 量 均 呈 微弱 的 增加 趋势 ， 
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最 “上 四 分 位 数 中 位 线 均值 下 四 分 位 数 最 
— 积 雪 覆盖 率 ZASE k A > 
值 值 


4800 F (a) 琼 乌 散 库 什 河流 域 ]45 4500 F (b) 木 扎 提 河 流域 70 
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图 5 6 一 9 月 雪线 高 度 与 平均 积 雪 和 覆盖 率 变化 
Fig. 5 Variation of snowline altitude and snow cover percentage from June to September 
81°10'E 81°40'E 
阿 克 牙 孜 河流 域 
Z 
3 | 
E 
Z 
S Le 
S 
81°10'E 81°40'E 

图 例 一 6H(3775m) — 8 月 (3936m) — 104 (3458 m) 雪线 高 度 /Im M 


一 7 月 (3858m) 一 9 月 (3692 m) 一 融雪 末期 雪线 (3967 m) 6332 1748 
图 6 瞬时 雪线 高 度 空间 分 布 (以 阿 克 牙 孜 河流 域 为 例 ) 


Fig.6 Spatial distribution of instantaneous snowline altitude (Take the Akyazi River basin for example) 


木 扎 提 河 流域 `. 阿 元 牙 孜 河流 域 和 玛 纳 斯 河流 域 年 温 (1.9 C) 最 低 , 玛 纳 斯 河流 域 多 年 年 平均 气温 
降水 量 和 夏季 降水 量 均 呈 减 少 趋势 ;各 流域 年 均 气 (3.3 %) 和 夏季 平均 气温 (8.7 “Cries BRS BE 
温和 夏季 平均 气温 均 呈 上 升 趋势 (图 8)。 其 中 , 木 。 河流 域 多 年 年 平均 降水 量 (485.8 mm) 和 夏季 平均 降 
扎 提 河流 域 多 年 年 平均 气温 (-8.9 色 ) 和 夏季 平均 气 。” 水 量 (306.3 mm) 最 少 , 阿 克 牙 孜 河流 域 多 年 年 平均 
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4700 + R2—0.1418 4400 上 R?-0.1871 
Š 0-94 E 0-59 
ag 4500- ,° we sc NE Ap 4200 上 -rr 
T 4300 L WT 4900 上 
4100 1 1 1 1 1 1 3800 1 1 1 1 1 Loa 
1990 1995 2000 2005 2010 2015 2020 1990 1995 2000 2005 2010 2015 2020 
年 份 年 份 
4400 F (c) 阿 克 牙 孜 河流 域 ui ume 4600 F (a) 玛 纳 斯 河流 域 
4200 上 —85 i y=6.3791x-8603.6 
E TS g 749 R°=0.3733 . 
iM 4000 上 VUE a Loans EN T | sueco 
wo |o oc ai 4200 - ee ae 
a 3800 |- x eee 
b RTT 
3600 上 4000 上 
3400 3800 
1990 1995 2000 2005 2010 2015 2020 1990 1995 2000 2005 2010 2015 2020 
年 份 年 份 
注 : 由 于 云层 遮挡 和 影像 缺失 等 原因 ,有 些 年 份 获取 不 到 融雪 末期 雪线 高 度 。 
图 7 融雪 末期 雪线 高 度 年 际 变化 
Fig.7 Interannual variation of snowline altitude at the end of melting season 
表 3 融雪 末期 雪线 高 度 与 气温 、 降 水 量 之 间 的 相关 性 
Tab.3 Correlation between snowline altitude at the end of melting season, temperature and precipitation 
夏季 气温 相关 性 年 气温 相关 性 夏季 降水 量 相关 性 年 降水 量 相关 性 
琼 乌 散 库 什 河流 域 0.485 0.352 -0.105 -0221 
木 扎 提 河流 域 0.558" 0.451 -0.232 -0.133 
阿 克 牙 孜 河流 域 0.503” 0.361 -0.136 -0.321 
玛 纳 斯 河流 域 0.391 0.366 -0.322 -0.371 


注 :*、** 分 别 表示 通过 0.05 和 0.01 的 显著 性 水 平 检验 。 


降水 量 (930.6 mm) 和 夏季 平均 降水 量 (596.1 mm ) 最 
多 。 人 研究 区 融雪 末期 雪线 高 度 与 气温 之 间 呈 正 相 
关 关 系 ,特别 是 与 夏季 气温 之 间 呈 显著 正 相 关 关 系 
(P< 0.05) ; 而 融雪 末期 雪线 高 度 与 降水 量 之 间 呈 较 
弱 的 负 相关 关系 ,这 表明 夏季 气温 是 影响 融雪 末期 
雪线 高 度 变化 的 主要 因素 。 
为 探究 研究 区 融雪 末期 雪线 高 度 对 气候 变化 
的 敏感 性 ,采用 多 元 线性 回归 方法 建立 了 融雪 末期 
雪线 高 度 与 夏季 气温 和 年 降水 量 回 归 模 型 ; 
SLA =À * aT -bP (7) 
式 中 :SL4 为 融雪 末期 雪线 高 度 ;和 A 为 常数 ;T 为 夏季 
气温 (C);P 为 年 降水 量 (mm);a.5 为 回归 系数 。4 
个 流域 回归 系数 如 表 4 所 示 ,以 阿 克 牙 孜 河 流域 为 
[n] , 当 夏 季 和 气温 升 高 或 降低 1 % 时 , 阿 克 牙 孜 河流 域 
的 融雪 末期 雪线 高 度 将 升 高 或 降低 76.2 m; 当年 降 


水 量 增加 或 减少 100 mm 时 ,该 流域 的 融雪 末期 雪线 
高 度 将 降低 或 升 高 28.7 m。 说 明 在 阿 克 牙 孜 河流 
域 ,夏季 气温 升 高 1 CC 引起 的 融雪 末期 雪线 高 度 上 
升 需要 增加 265 mm 的 降水 量 才 能 弥补 。 琼 乌 散 库 
什 河 流域 夏季 气温 回归 系数 最 大 (82.9) ,意味 着 该 
流域 融雪 末期 雪线 高 度 变化 对 夏季 气温 响应 最 为 
敏感 。 玛 纳 斯 河流 域 年 降水 回归 系数 最 大 (0.38)， 
说 明 该 流域 融雪 末期 雪线 高 度 变 化 对 年 降水 量 响 
应 最 为 敏感 。 


4 讨论 


4.1 积 雪 制 图 与 雪线 高 度 提取 的 不 确定 性 
积 雪 制图 和 雪线 高 度 提取 的 不 确定 性 主要 来 
自 以 下 4 个 方面 :(1) 利用 光学 遥感 影像 提取 山区 
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图 8 融雪 末期 雪线 高 度 与 气温 和 降水 量 的 年 际 变化 


Fig.8 Interannual variations of temperature and prec 


表 4 融雪 末期 


ipitation between snowline altitude at the end of melting season 


雪线 高 度 回 归 模型 系数 


Tab.4 Regression model coefficients of snowline altitude at the end of melting season 


- 2 一 = 相关 系数 (R) 显著 性 水 平 
琼 乌 散 库 什 河流 域 4220 82.9 0.355 0.515 P< 0.05 
木 扎 提 河 流域 4294 50.7 0.238 0.559 P< 0.01 
阿 克 牙 孜 河流 域 4040 76.2 0.287 0.548 P< 0.01 
玛 纳 斯 河流 域 4009 62.4 0.364 0.520 P« 0.05 


积 雪 信息 时 ,地 形 起 伏 产 生 的 山体 阴影 会 影响 积 雪 
范围 提取 的 准确 性 ,这 种 影响 也 被 称 为 地 形 效应 。 
通过 地 形 校正 * 可 以 较 好 地 消除 地 形 效应 对 山区 
积 雪 识别 的 影响 。 此 外 ,如 果 是 像 元 尺度 雪线 高 度 
的 提取 ,还 应 定量 分 析 坡 向 对 雪线 高 度 的 影响 。 本 
研究 专注 的 是 区 域 性 雪线 高 度 的 提取 。 该 方法 求 
解 得 到 的 唯一 高 程 值 代 表 了 一 个 区 域内 的 综合 雪 
线 高 度 。 因 而 ,本 研究 没有 考虑 地 形 效 应 以 及 坡 向 
等 微 地 形 因子 的 影响 。(2) 云 和 云 阴 影 的 存在 也 是 
影响 光学 遥感 中 积 雪 信息 提取 的 主要 误差 来 源 。 
本 文 基于 Landsat 系列 数据 ,采用 瞬时 雪线 去 云 和 时 
序 融 合 的 方法 进行 去 云 处 理 。 而 Sentinel-1 搭载 了 


合成 孔径 雷达 (SAR ) ,该 数据 具有 全 天 时 、 全 天 候 和 
穿 透 性 的 特点 ,不 受 天 气 和 太阳 光照 条 件 的 影响 ; 
Sentinel-2 数 据 为 多 光谱 卫星 成 像 ,该 数据 具有 较 高 
的 时 间 分 辨 率 (5 d) 。 在 后 续 的 进一步 研究 中 可 以 
结合 Sentinel-1 和 Sentinel-2 数 据 ,增强 云 下 积 雪 识 
别 的 能 力 。(3) 对 于 和 森林 覆盖 区 积 雪 的 提取 精度 较 
差 , 树 木 的 冠 层 会 挡住 地 面部 分 的 积 雪 信息 。 基 于 
NDSI NDFSI (Normalized Difference Forest Snow In- 
dex) , NDVI(Normalized Difference Vegetation Index) 
的 多 指标 积 雪 识 别 方法 ””, 可 有 效 提取 林 下 积 雪 
言 息 。 但 融雪 末期 雪线 高 度 要 远 高 于 森林 植被 区 ， 
所 以 本 人 研究 没有 考虑 森林 的 影响 ,(4) 因 大 部 分 裸 
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露 冰 川 表 面 不 清洁 ( 表 厂 覆盖 ) 反 照 率 明显 低 于 积 
雪 , 因 而 不 清洁 的 裸露 冰川 和 非 季 节 性 积 雪 可 以 通 
过 光学 遥感 (NDSI 阅 值 ) 进 行 区 分 。 然 而 ,融雪 末期 
雪线 高 度 的 提取 ,也 应 该 注重 区 分 非 季 节 性 积 雪 与 
表面 清洁 的 冰川 ;例如 可 利用 两 者 在 近 红 外 波段 上 
所 表现 出 的 不 同 特 征 ,根据 近 红 外 波段 影像 直方 图 
分 类 出 冰川 和 积 雪 ””。 特 别 是 随 着 Sentinel-1 等 
星 载 合成 孔径 雷达 的 出 现 ,使 冰川 和 积 雪 的 区 分 更 
具 可 能 性 。 

本 研究 中 的 融雪 末期 雪线 高 度 提 取 方 法 仅 适 
合 于 区 域 (或 者 流域 ) 尺 度 的 雪线 高 度 提 取 。 但 若 
注重 分 析 像 元 尺度 雪线 高 度 的 时 空 差 异性 (比如 坡 
向 .坡度 等 微 地 形 因 子 的 影响 ) ,应 采用 像 元 尺度 雪 
Zea BERET EP 
42 融雪 末期 雪线 高 度 的 空间 分 异 

本 研究 发 现 , 近 30 a 来 研究 区 域 融雪 末期 雪线 
高 度 呈 上 升 趋势 ,这 与 天 山 物 质 平 衡 处 于 负 平衡 状 
态 的 研究 结论 是 一 致 的 后 。 从 研究 区 域 每 年 的 融 
雪 末 期 雪线 高 度 来 看 , 琼 乌 散 库 什 河流 域 总 体 融 雪 
末期 雪线 高 度 明 显 高 于 其 他 3 个 流域 , 阿 克 牙 孜 河 
流域 融雪 末期 雪线 高 度 处 于 最 低 水 平 。 造 成 这 种 
雪线 高 度 分 布 特征 的 主要 原因 为 : 琼 乌 散 库 什 河流 
域 位 于 天 山南 坡 ,天 山南 坡 为 阳 坡 ,受到 太阳 辐射 
量 较 多 ,造成 气温 相对 于 北 坡 较 高 的 ; 且 该 流域 地 
势 总 体 比 较 高 , 受 山 体 海拔 影响 ,所 以 该 地 区 融雪 
末期 雪线 高 度 最 高 ; 阿 克 牙 孜 河 流域 位 于 天 山北 
坡 , 受 大 西洋 和 北冰洋 的 水 汽 影响 中 ,降水 量 相 对 
充沛 ,太阳 辐射 量 少 ,气温 相对 较 低 ,因而 阿 克 牙 孜 
河流 域 融 雪 末 期 雪线 高 度 最 低 。 
4.3 气候 变化 对 融雪 末期 雪线 高 度 的 影响 

雪线 作为 一 种 气候 标志 线 , 其 高 度 受气 候 因素 
的 影响 同时 也 会 反映 气候 的 变化 。 现 有 的 研究 所 所 
表明 :在 全 球 气候 变 暖 的 影响 下 冰川 面积 不 断 奢 
缩 ,雪线 高 度 持 续 升 高 ,这 与 本 文 的 研究 相 一 致 。 
本 研究 发 现 , 夏 季 气 温 是 影响 融雪 末期 雪线 高 度 变 
化 的 主要 因素 ,以 往 有 关 高 山区 融雪 末期 雪线 或 者 
冰川 平衡 线 变 化 的 研究 也 得 出 类 似 结论 :Tang 等 ” 
基于 MODIS 数据 从 30 km 格 网 尺度 也 发 现 了 
2001 一 2016 年 亚洲 高 山区 融雪 末期 雪线 高 度 总 体 
量 上 升 趋势 ,特别 是 在 天 山地 区 和 喜马拉雅 地 区 上 
升 明显 ,夏季 气温 是 主要 影响 因素 ; 陈 安 安 等 E 
宁 练 等 所 :发现 ,夏季 气温 是 影响 冰川 平衡 线 变 化 的 


主导 因子 。 雪 线 高 度 与 气温 呈正 相关 关系 ,气温 升 
高 会 促进 雪线 高 度 的 升 高 ;雪线 高 度 与 降水 量 呈 负 
相关 关系 ,降水 量 增加 会 引起 雪线 高 度 的 降低 。 在 
后 续 的 研究 中 可 以 基于 遥感 提取 的 长 时 序 雪线 高 
度 .气温 和 降水 等 气象 数据 ,利用 机 噩 学 习 方 法 构 
建 雪 线 高 度 变化 的 预测 模型 。 


5 结论 


AS SEF GEE Bit OF S All Landsat 影像 ,发 
展 了 融雪 末期 区 域 雪 线 高 度 遥 感 提 取 模 型 ,提取 了 
天 山 4 个 流域 近 30 a 长 时 间 序 列 的 融雪 末期 雪线 高 
度 ;并 对 4 个 流域 融雪 末期 雪线 高 度 的 时 空 变化 特 
征 及 其 与 气温 、 降 水 之 间 的 关系 进行 分 析 , 得 出 以 
下 结论 : 

CL) 模型 提取 的 融雪 末期 雪线 高 度 与 Sentinel- 
2 提取 的 融雪 末期 “最 小 化 ” 积 雪 范围 具有 较 高 的 一 
致 性 ,平均 总 体 精度 为 91.6%;Kappa 系数 平均 值 为 
0.9。 利 用 该 提取 模型 可 准确 提取 流域 内 融雪 末期 
区 域 雪 线 高 度 , 可 有 效应 用 于 流域 尺度 .长 时 间 序 
列 雪线 高 度 的 遥感 监测 。 

(2) 研究 区 1991—2021 年 融雪 末期 雪线 高 度 
年 际 波动 较 大 ,但 总 体 上 呈现 出 上 升 趋势 。 阿 克 牙 
孜 河流 域 融 雪 末 期 雪线 高 度 上 升 速率 最 慢 (2.8 me 
a^) , 玛 纳 斯 河流 域 融 雪 末 期 雪线 高 度 上 升 速率 最 
快 (6.4 mea’) 

(3) 研究 区 融雪 末期 雪线 高 度 主要 受气 温 的 影 
响 , 与 夏季 气温 之 间 具 有 显著 的 正 相关 关系 ,与 降水 
量 之 间 具 有 较 弱 的 负 相 关 关 系 。 当 夏季 气温 升 高 
(降低 )1 % 时 ,流域 融雪 末期 雪线 高 度 约 升 高 (降低 ) 
50.7-82.9 m; 当年 降水 量 增加 (减少 )100 mm 时 , 流 
域 融 雪 末期 雪线 高 度 将 约 降 低 ( 升 高 )23.8~36.4 mo 

(4) GEE 平台 强大 的 存储 与 计算 能 力 等 ,大 幅 
提高 了 雪线 高 度 遥 感 提取 的 效率 ,为 长 时 序 .大 范 
于 遥感 雪线 高 度 监测 提供 了 强 有 力 的 计算 资源 
保障 。 
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Monitoring of snowline altitude at the end of melting season in Tianshan 
Mountains from 1991 to 2021 


WANG Jingwen', TANG Zhiguang, DENG Gang, HU Guojie, SANG Guoging' 
(1. Hunan Provincial Key Laboratory of Geo-Information Engineering in Surveying, Mapping and Remote Sensing, 
Hunan University of Science and Technology, Xiangtan 411201, Hunan, China; 2. Northwest Institute of 
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Abstract: The study of snowline altitude changes at the end of the melting season is helpful to predict future 
trends of snow systems in order to understand the regional and global climate changes. The remote sensing 
extraction method of region snowline altitude was developed based on the Google Earth Engine and Landsat 
satellite data. From 1991 to 2021, the snowline altitude at the end of the melting season in the four basins of the 
Tianshan Mountains was extracted. The variation characteristics of snowline altitude and its relationship with 
meteorological factors were carefully and accurately estimated. The results were as follows: (1) the extended 
snowline altitude at the end of the melting season correlated well with the snow cover extent (minimum) at end of 
the melting season, which was extracted by Sentinel-2. The overall accuracy was 91.6%, and the Kappa coefficient 
was higher than 0.9. The regional snowline altitude at the end of the melting season was accurately obtained with 
this model. (2) Generally, the snowline altitude at the end of the melting season in the study area showed an 
obvious increasing trend over the last 30 years. The increase rate was between 2.7 m-a' and 6.4 m-a''. The rise rate 
of the snowline altitude in the Manas River Basin was the fastest (6.4 m + a`’), while the rise rate of the snowline 
altitude in the Akeyazi River Basin was the slowest (2.7 m-a’'). (3) The summer temperature was the main factor 
that affected the change in the snowline altitude at the end of the melting season in the study area (P < 0.05), 
while the effect of precipitation was relatively weak. 


Keywords: snowline altitude; Landsat; remote sensing monitoring; climate change; Tianshan Mountains 


